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Abstract 

The Dirac equation for an electron in two spatial dimensions in the Coulomb and ho- 
mogeneous magnetic fields is an example of the so-called quasi-exactly solvable models. 
The solvable parts of its spectrum was previously solved from the recursion relations. 
In this work we present a purely algebraic solution based on the Bethe ansatz equa- 
tions. It is realised that, unlike the corresponding problems in the Schrodinger and the 
Klein-Gordon case, here the unknown parameters to be solved for in the Bethe ansatz 
equations include not only the roots of wave function assumed, but also a parameter 
from the relevant operator. We also show that the quasi-exactly solvable differential 
equation does not belong to the classes based on the algebra sfe- 
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1. Recently a new type of spectral problem, so-called quasi-exactly solvable model 
(QESM), was discovered by physicists and mathematicians [PJ-|§]. This is a special class 
of quantum-mechanical problems for which analytical solutions are possible only for parts of 
the energy spectra and for particular values of the fundamental parameters. The reason for 
such quasi-exactly solvability is usually the existence of a hidden Lie-algebraic structure [@- 
@]. More precisely, quasi-exactly solvable (QES) Hamiltonian can be reduced to a quadratic 
combination of the generators of a Lie group with finite-dimensional representations. 

The first physical example of QESM in atomic physics is the system of two electrons 
moving in an external oscillator potential discussed in || [TIJ. The authors of these works 
apparently were unaware of the mathematical development in QESM. Later, several physi- 
cal QESM were discovered, which include the two-dimensional Schrodinger [jUL the Klein- 



Gordon |fL2| , and the Dirac equations of an electron moving in an attractive/repulsive 
Coulomb field and a homogeneous magnetic field. The essential features shared by all these 
above examples are as follows. The differential equations are solved according to the stan- 
dard procedure. After separating out the asymptotic behaviors of the system, one obtains 
an equation for the part which can be expanded as a power series of the basic variable. But 
instead of the two-step recursion relations for the coefficients of power series so often encoun- 
tered in exactly solvable problems, one gets three-step recursion relations. The complexity 
of the recursion relations does not allow one to determine the energy spectrum exactly from 
the normalisability of the eigenfunctions. However, one can impose a sufficient condition for 
normalisability by terminating the series at a certain order of power of the variable; i.e. by 
choosing a polynomial. By doing so one could obtain exact solutions to the original problem, 
but only for certain energies and for specific values of the parameters of the problem. These 
parameters, namely, are the frequency of the oscillator potential and the external magnetic 
fields. 



In [[TJJ a systematic and unified algebraic treatment was given to the above-mentioned 
systems, with the exception of the Dirac case. This was made possible by realising that 
these systems are governed essentially by the same basic equation, which is quasi-exactly 
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solvable owing to the existence of a hidden sl 2 algebraic structure. This algebraic structure 
was first realised by Turbiner for the case of two electrons in an oscillator potential [HJ . In 
this algebraic approach, analytic expressions of the solvable parts of the energy spectrum 
and the allowed parameters were expressible in terms of the roots of a set of Bethe ansatz 
equations. 



In this paper we would like to extend the method of [14] to the planar Dirac equation 
of an electron in the Coulomb and magnetic fields. It turns out the a set of Bethe ansatz 
equation can also be set up in this case. However, unlike the systems considered in [H, 
here the unknown variables in the Bethe ansatz equations involved not only the roots of 
the wave functions assumed, but also a parameter from the relevant operator. We also 
demonstrate that the Bethe ansatz approach yields the same spectrum as that obtained by 
solving recursion relations. Finally, we show that the quasi-exactly solvability of this system 
is not related to the sl 2 algebra. 

2. In 2+1 dimensions the Dirac algebra 

h", y} = 2<r , <r = dia g (i, -1, -1) (i) 

may be represented in terms of the Pauli matrices as 7 = 03, 7 fc = iaj., or equivalently, 
the matrices (ai,a 2 ) = 7°(7 1 ,7 2 ) = (— cr 2 ,cr 1 ) an d /3 = 7 . Then the Dirac equation for 
an electron minimally coupled to an external electromagnetic field has the form (we set 

c = h = i) 

(idt-H D )V(t,r) = 0, (2) 

where 

H D = aP + (3m - eA° = a x P 2 - o 2 P x + a 3 m - eA° (3) 

is the Dirac Hamiltonian, = —idk + eAk is the operator of generalized momentum of 
the electron, the vector potential of the external electromagnetic field, m the rest mass 
of the electron, and — e (e > 0) is its electric charge. The Dirac wave function ty(t,r) 



is a two-component function. In an external Coulomb field and a constant homogeneous 
magnetic field B > along the z direction, the potential assume the following forms in 
the symmetric gauge 

A°(r) = Ze/r (e > 0), A x = -By/2, A y = Bx/2 . (4) 

We assume the wave functions to have the form 

tt(t, x) = exp(-iEt)i/}i(r, <p) , (5) 



where E is the energy of the electron, and 

F(r)e 



(6) 



with integral number I. The function ipi(r,(p) is an eigenfunction of the conserved total 
angular momentum J z = L z + S z = —id/dip + cr 3 /2 with eigenvalue j — I + 1/2. It should 
be reminded that I is not a good quantum number. Only the eigenvalues j of the conserved 
total angular momentum J z are physically meaningful. 

By putting Eq.(|5|) and (|6]) into (0), and taking into account of the equations 

,. ±1 p.-^(» ± (1£-*)). p, 

we obtain 

dF + \ eBr\ „ / Za\ ^ 

*"(^ + — )F+(s + m + -)G = 0, (8) 

where a = e 2 = 1/137 is the fine structure constant. In a strong magnetic field the asymp- 
totic solutions of F(r) and G(r) have the forms exp(— eBr 2 /4) at large r, and r 7 with 



J — y(l + 1/2) 2 — (Za) 2 for small r. One must have < 1/2, otherwise the wave function 
will oscillate as r — > when I = and I = —1. 
Let us assume 

F{r) = r 7 exp(-eBr 2 /4) Q(r), G(r) = r 7 exp(-eBr 2 /4) P{r) . (10) 



In |13| we showed that parts of the spectrum could be analytically solved for by imposing 
the sufficient condition that Q(r) and P(r) be polynomials, thus showing that the system 
belongs to the QESM. The spectrum was solved in [jn| from the recursion relations for the 
coefficients in the series expansion in Q and P. In this paper, we will show that the same 
spectrum can also be obtained in a purely algebraic way This is achieved by the method of 
factorisation which leads to a set of Bethe ansatz equations |13|, [14]] . 

Substituting Eq.flTOD into Eq.(§) and (^j) and eliminating P{r) from the coupled equations, 
we have 



dr 2 



+ 



r E + m + Zat/r 



^ E 1 
dr 



m 



lEZa I + \ 



+ 



7 



Zajr 



E + m + Zajr 



eisr 



Z + 1/2' 



Q(r) = 



where T = Z + 1/2 + 7. Once Q(r) is solved, the form of P(r) is obtainable from Eqs.(§) and 
(jig). If we let x = r/l B ,l B = 1/VeB, Eq.flll]) becomes 



dx 2 



+ 



2 7 Za 

— — x ~\~ 

x x((E + m)l B x + Za) 



+ 



X 



x^ 



Za(Z + 1/2 -7) 



x 2 [(E + m)l B x + Za] (E + m)lBX + Za 
Eq.([E|) can be rewritten as 



<ir 2 



2/3 



x 



X 



X + Xq 



— + e + - 

dec oc 



■(T + l) 
Q(x) = . 

g(x) = 



(12) 



(13) 



Here = 7 + 1/2, x = Za/[(E + m)l B ), e = (E 2 - m 2 )l 2 B - (T + 1), b = b + L/x , 
b = 2EZal B , L = (I + 1/2 — 7), and c = x + L/x . On expressing l B in the expression of 
e in terms of x , we get 



E — m f Za 



r + r 



(14) 



E + m \ Xq J 

It is obvious that the energy E is determined once we know the values of e and Xq. The 
corresponding value of the magnetic field B is then obtainable from the expression l B = 
Za/[(E + m)xo\- Solution of xq is achieved below by means of the Bethe ansatz equations. 



3. We observe that the problem of finding the spectrum for Eq.flTB"!) is equivalent to 
determining the eigenvalues of the operator 

H=-^--(y.-x--±-)±-*+-?-. (15) 
dx 2 \ x x + xq ) dx x x + xq 



We want to factorise the operator dig) in the form 

H = a + a + e. (16) 

The eigenfunctions of the operator H at e = must satisfy the equation 

a Q(x) = . (17) 

Suppose polynomial solution exist for Eq.flTSp, say Q equals a non- vanishing constant, or 
Q — IYk=i( x ~ x k), where x^ are the zeros of Q, and n is the degree of Q. In the case where 
Q is a constant (which may be viewed as corresponding to n = 0), the operators a and a + 
have the form 

d , d [2/3 1 \ 

a = — , a = — x 



d 


2(3 


dx 


1 X 



dx dx \ X X + Xq 

If Q — IYk=i( x ~ x k), a, and a + will assume the form 

d A 1 



and 



(XX fc— 1 ^ 



d [2(3 1 \ ™ 1 



— x — 



E z—r ■ (2°) 



dx \ 

We now substitute the forms of a and a + into Eq.(J16|) and compare the result with 
Eq.(^). This leads to conditions that must be satisfied by the various parameters and the 
roots Xfc's. For constant Q (n = 0), one has 

e = b = c = 0. (21) 

The fact that c = implies 

4 = -L. (22) 
6 



For n > 1, one gets 

L n 1 

6 + - = 2/3 ]T-, e = n, (23) 

•£0 fc=i 

L n 1 

^o + - = ' (24) 

2/3 1 ™ 1 , 

— - — 2^ = 0, fc = l,...,n. (25) 

Eqs. (|22"D , fl24 ) and (|25 ) constitute the set of n + 1 Bethe ansatz equations relevant to this 



Dirac system, which involve n + 1 unknown parameters {xo,Xi, . . . ,x n } . It is worthwhile 
to note that, unlike the corresponding equations in the Schrodinger and the Klein-Gordon 
case discussed in [|4j], this set of Bethe ansatz equations involved not only the roots x^s, 



but also a parameter x from H. From the second equation in Eq. fl23|) we get 

E 2 -m 2 = ^iT + n+l) . (26) 

Since — l/2<r<0for Za < 1/2 0], we see from Eq. (p6|) that the solvable parts of the 
spectrum must satisfy \E\ > m. 

So we see that the solution of the solvable parts of the spectrum E boils down to solving 
the Bethe ansatz equations for xq in the differential operator, and the roots Xk {k = 1, . . . ,n) 
of Q(x). Once the value of xq for each order n = e is known, the energy E is given by 
Eq.(|n]). The corresponding magnetic field B is then determined from the definition of 60, or 
from Eq . (|26|) . The Bethe ansatz equations thus provides a systematic solutions of the QES 
spectrum. Of course, as the order of the degree of Q increases, analytical solutions of the 
Bethe ansatz equations becomes difficult, and one must resort to numerical methods. 

4. In what follows we shall show the consistency of the solutions by the Bethe ansatz 



approach and that by the recursion relations presented in jT3| for the first three lowest orders 
(n = 0,1,2) in Q. Instead of solving for xq, our strategy is to eliminate it in Eq.(TTJ]) by 
means of the equations (p2h-Q25|) so as to obtain an equation obeyed by E for each order 
of Q. This equation is then compared with the corresponding equation obtained from the 
recursion relations as presented in Jl3| . 



From Eq.([H|) and (|22| ) we have Xq = —L and e = when Q is a constant. Substitute 
these values of Xq and e into Eg. (|T4"D , and using the fact that TL = (Za) 2 , we obtain the 
corresponding value of E as 



E 



rn 



2(1 + j + 1) ' 



(27) 



This is the result presented in [13]. The corresponding allowed value of the magnetic field 
B is then obtained from Eq.(|2BJ) and (P7|). 



For n = 1, we find from Eqs.(|l4|), ( [23] 



r + 2 



b + — = 

x 



- Xi 



1 and (II) that 

E-m (Za) 2 
E + m Xq 
2£ 

i 

Xi 

L 

Xq H , 

Xq 



1 







Eq.( 



and 



Xi Xi + Xq 

imply X\ — &o — Xq. Substituting x\ into Eq.( 

-i -l 



(28) 
(29) 
(30) 
(31) 



we obtain 



x 



o 



E + m 
2E(Za)' 



(32) 



Then from Eq.fl3"2"|) and Eq.i 



4(r + i)- 



we get 

r 

z 2 a 2 



E 2 + AEm + 



r 



(Za) 



-m 



(33) 



The energy E can be solved from Eq . (|3~3D by the standard formula, after which the magnetic 
field is determined from Eq. (p6|) . Eq.(|33|) does not resemble the one obtained from recursion 
relation in [jl3| . However, on multiplying Eq. ( |33| ) by V + 1 and making use of the fact that 
(Za) 2 = T(r — 27), we can show, after some algebra, that Eq.(^) is equivalent to the 



corresponding equation given in |I3 



Finally we consider the case for n = 2. We have Eq. fll4|) with e = 2, together with 
Eqs.(||), (||) and © in the forms 

E-m(Za) 2 



r + 3 



E + m Xq 
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(34) 



L 2(3 2(3 

On H = 1 



X %l %2 

1 1 L 
; + ; = x + — , 

X\ + X X 2 + X Xq 

2(3 1 2 n 
xx - ■ = , 

X\ x x + x x 2 — Xx 

^ 1 2 n 
x 2 - ■ = . 

x 2 x 2 + X Xx~ x 2 



(35) 
(36) 
(37) 
(38) 



From these equations we find xx + x 2 = b — xo and xxx 2 = 2(3xo(bo—xo)/(boXo + L). Putting 

these expressions into Eq. (|36|) and using the fact that r = 2(3 + L — 1, we arrive at 

b TL 



(b 2 - 2p) x 2 + b Tx + [b 2 (L — l) — L {2(3 + 1 



+ 







(39) 



Now multiplying Eq . (|3~9D by T, using TL = (Za) 2 , and expressing bo, Ib, and 1/xq in terms 
of E, we get finally 



4(2r + 3) - 



+ U2 



(Za 



(Za) 



+ 



(Za 



6r + 2( 7 + 1) + 

2(7 + 1)- 

6r + 2( 7 + 1) + 

2(7 + 1) 



(2 7 +i)r 
r + 3 

(2j + l)T 

r + 3 
(2 7 + i)r 



E 6 



E 2 m 



r + 3 
(27 + i)r 



r + 3 



Em 2 



m=0 



(40) 



Again, this equation does not look the same as that obtained from the recursion relations. 
But we can show they are in fact equivalent as they differ only by a multiplicative factor 

(r + i)(r + 2). 



5. We now demonstrate that the QES equation (|13|) cannot be represented as bilinear 
combination of the generators of the sl 2 algebra. The question of whether there exists non- 
s/2-based one-dimensional QESM was first posed in M in which all s/2-based QESM are 



classified. The first example of such a kind was given in ||16|| , which presents a potential 
arising in the context of the stability analysis around the kink solution for 4 -type field 
theory in 1 + 1 dimensions. 

To show that Eq . (|i~3"D is also not generated by the sl 2 algebra, let us rewrite it as 



d 2 



x 2 + x x) — — + \x^ + x x 2 + {l-2(3)x- 2(3x Q 
' ax 1 L 



d 

dx 



-ex 2 



+ (c- b - ex ) x - bx \ Q(x) = . (41) 



Turbiner has shown that all s/ 2 -based second order QES differential equations can be cast 
into the form 

rPO dO 

~ Pi{x) i^ + Ps{x) i^ + {P2[x) ~ A) Q = ' (42) 

where 

Pi(x) = a ++ x 4 + a +0 x 3 + + a 00 ) x 2 + a -X + a , 

P 3 (x) = 2(2j-l)a ++ x 3 +l(3j-l)a +0 + b + ]x 2 

+ [2j (a+_ + a o) + a 00 + b ] x + ja - + 6- , 
P 2 (x) = 2j (2j - 1) a ++ x 2 + 2j (ja +0 + b + ) x + a 00 f + b j . (43) 

Here a^'s and fo^'s (A;, / = +,0,—) are constants, and j is a non-negative integer or half- 
integer. Eq.([4"2D corresponds to the eigenvalue problem 



HQ = XQ , H = - £ a k iJ k J l + ^ > ( 44 ) 



fe,i=+,0,- fc=+,0,- 
k>l 



which has a polynomial solution of power 2j in x. Here J s are the generators of sl^- 

J + =x 2 ^--2jx, J° = x±-j, J~ = ^. (45) 

Comparing Eqs.(f4"l"D and ( ^2[ ) we find that the two equations are inconsistent with each 
other. For instance, the coefficient of x 4 in P 4 requires a ++ = 0, whereas the coefficient of 
x 3 in P3 implies 2(2j — l)a ++ = 1, which gives a non-vanishing a ++ for positive integral and 
half-integral values of j. This shows that Eq.([13|) is not s/ 2 -based. 

6. In conclusion, we have given an algebraic solution to the planar Dirac equation of 
an electron in the Coulomb and magnetic fields. The relevant Bethe ansatz equations are 
presented. Unlike the corresponding equations in the Schrodinger and the Klein-Gordon case 
discussed in jTJJ] , the unknown variables in this set of Bethe ansatz equations include not only 
the roots of the polynomial assumed, but also a parameter from the QES differential operator. 

10 



Equivalence between this approach and that by the recursion relations is demonstrated. 
Finally, we show that the QES equation for this problem does not belong to any of the 
classes based on the sfa algebra. 
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